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ABSTRACT 

Hierarchical galaxy mergers will lead to the formation of binary and, in the case of a subsequent merger 
before a binary coalesce, triple supermassive black holes (SMBHs), given that most massive galaxies harbor 
SMBHs. A triple of SMBHs becomes visible as a triple active galactic nucleus (AGN) when the BHs ac- 
crete large amounts of gas at the same time. Here we report the discovery of a kpc-scale triple AGN, SDSS 
J1027+1749 at z = 0.066, from our systematic search for hierarchical mergers of AGNs. The galaxy contains 
three emission-line nuclei, two of which are offset by 450 and 110 km s" 1 in velocity and by 2.4 and 3.0 
kpc in projected separation from the central nucleus. All three nuclei are classified as obscured AGNs based 
on optical diagnostic emission line ratios, with black hole mass estimates M, > 10 8 M Q from stellar velocity 
dispersions measured in the associated stellar components. Based on dynamical friction timescale estimates, 
the three stellar components in SDSS J1027+1749 will merge in ~ 40 Myr, and their associated SMBHs may 
evolve into a gravitationally interacting triple system in < 200 Myr. Our result sets a lower limit of ~ 5 x 10~ 5 
for the fraction of kpc-scale triples in optically selected AGNs at z ~ 0.1. 

Subject headings: black hole physics — galaxies: active — galaxies: interactions — galaxies: nuclei — 
galaxies: stellar content 



1. INTRODUCTION 

G alaxies are thought to be b uilt up hierarchically via merg- 
ers dToomre & Toomrell 19721) . Because most massive galax- 
ies are b elieved to harbor a central sup ermassive black hole 
(SMBH; [Kormendy & Richstondll995f) . galaxy mergers will 
result in the formation of binary SMBHs through dynam- 
ical f ricti on dBeg elman et al. 1980; Milosavljevic & Merritt 
1200 It lYul 120021) . When a binary does not coalesce before 
a subsequent merger with a third galaxy, a system of three 
gravitationally interacting SMBHs is expected to form (e.g., 
IValtonenll 19961) . Triple SMBHs may hold important clues to 
our general understanding of massive galaxy formation. For 
example, numerical simulations suggest that triple SMBHs 
scour out cores in stellar bulges with mass deficits one or 
two times of the total mass of the BHs and sizes larger than 
those formed around binary SMBHs. This process may be 
responsible for the large core s observed in some ma ssive el- 
liptical galaxies such as M87 dHoffman & Loebl2007l) . Triple 
SMBHs also provide a unique astrophysical laboratory to 
study the chaotic dyna mics of three-body interactions i n Gen- 
eral Relativity (e.g., iBlaes et al.ll2002t iMerrittl 12006b . Nu- 
merical simulations of the dynamics of triples suggest that 
their encounter may lead to a merger of all three BHs, for- 
matio n of a highly eccentric binary, or ejection of th ree free 
BHs dlwasawaet al. 2006; lLousto & Zlochowej2008ft . Triple 
SMBHs exhibit phases of very high eccentricity in the in- 
ner binary, producing intense bursts of gravitational radiation 
which will be within the sensitivity range of forthcoming pul- 
sar tim ing arrays and the Laser I nterferometer Space Antenna 
(e.g JAmaro-Seoane et alll2010l) . 

Despite their significant scientific merit and intense theoret- 

1 This paper includes data gathered with the Apache Point Observatory 
3.5-meter telescope, which is owned and operated by the Astrophysical Re- 
search Consortium. 

2 Harvard-Smithsonian Center for Astrophysics, 60 Garden St., Cam- 
bridge, MA 02138 

3 Einstein Fellow 

4 Princeton University Observatory, Princeton, NJ 08544 



ical interest, direct observational evidence for gravitationally 
interacting triple SMBHs is still lacking, however, because 
their typical separation (less than a few parsecs) is too small 
to be resolved at cosmological distances. In this Letter, we 
report the probable discovery of a kpc-scale triple of SMBHs, 
at a separation large enough that the components can be eas- 
ily resolved using current facilities, yet small enough for the 
system to be dynamically interesting. Using a simple stel- 
lar dynamical friction argument, we estimate that the system 
will form a bound SMBH triple in < 200 Myr. A triple of 
SMBHs becomes visible as a triple active galactic nucleus 
(AGN) when the BHs accrete large amounts of gas at the same 
time - a proc ess which is thou ght to be common in gas-rich 
mergers (e.g., Hernquist 1989). While kpc-scale triple AGNs 
are supposed to be present in recurrent galaxy mergers, only 
one candidate is known 5 . It is a triplet of emission-line nu- 
clei (with two offset nuclei at 5.1 and 8.4 kpc in projection 
from the primary nucleus) in the minor merger NGC 334 1 at 
z ~ 0.03 serendipitously discovered by Bart h et al.l (120081) . 

To quantify the frequency of kpc-scale binary and triple sys- 
tems of SMBHs, we are conducting a systematic search for 
hierar chical mergers of AGN s in the Seventh Data Release 
(DR7; lAbazaiian et alj|2009l) of the Sloan Digital Sky Survey 
(SDSS: lYorketal.ll2000l) . We have selected a sample of 1286 
AGN pairs with projected separations r p < 100 kpc and line- 
of-sight velocity offsets Av < 600 km s" 1 dLiu et alJl201 lbl) . 
Ninety-two AGN pairs in the sample have r p < 10 kpc and 
z < 0. 16. Seven of these 92 pairs have a third galaxy or nuclei 
with r p < 10 kpc from both of the first two nuclei and r-band 
magnitude differences < 4 mag. Because of its proximity 
with the first two nuclei, none of the third nuclei has an SDSS 
spectrum due to the finite size of SDSS fibers 6 . To determine 

5 The re is one probable phys ical triple quasar known, QQ 1429-008 at 
z ~ 2.1 iDjorgovski et al. 2007), but the projected separations between the 
quasar components are much larger, ~ 30^50 kpc. 

6 The first two nuclei both have an SDSS spectrum because they were 
observed on overlapping spectroscopic plates. 
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FIG. 1.— Left: SDSS gri-color composite image of SDSS J1027+1749, 
shown with a 1' X 1' F OV. Contours indic ate radio flux densities from the 
FIRST 1.4 GHz survey iBecker et al. 1995). The galaxy has three nuclei A, 
B, and C, and a strong tidal feature to the northeast of B. Nucleus B (C) has 
a projected separation of 3.0 (2.4) kpc from A. Right: SDSS r-band image of 
the central 30" X 30" region. 

whether the third nucleus is also active and physically close 
to the first two nuclei, we are conducting spatially resolved 
spectroscopy for these kpc-scale triple AGN (hereafter, triple 
AGN, for short) candidates. Here, we report our initial results 
on the discovery of a triple AGN in a hierarchical merging 
system, SDSS J1027+1749. Shown in Figure Q] the galaxy 
contains three nuclei as seen in its SDSS image, with B and C 
at projected separations of 3.0 and 2.4 kpc from A. The nuclei 
A and B have SDSS spectra. We have conducted optical slit 
spectroscopy for all three nuclei. We describe our observa- 
tions and data analysis in ^2j and present physical measure- 
ments of the three nuclei in §|3] We discuss implications of 
our results and conclude in §|4] A ACDM cosmology with 
VL m = 0.3, f^A = 0.7, and h = 0.7 is assumed throughout. 

2. OBSERVATIONS AND DATA ANALYSIS 

We obtained slit spectra for SDSS J1027+1749 on the 
nights of 2011 March 2 and 6 UT using the Dual Imaging 
Spectrograph (DIS) on the Apache Point Observatory 3.5 m 
telescope. The sky was mostly non-photometric, with seeing 
ranging between 0."8 and 2."0 with a median around l."2. 
DIS has a 4' x6' FOV with a pixel size of 0."414. We adopted 
a l."5 x 6' slit with the B1200+R1200 gratings centered at 
510 and 700 nm. The spectral coverage was 450-560 (640- 
760) nm with an instrumental resolution of Cinst w 60 (30) 
km s" 1 and a dispersion of 0.62 (0.58) A pixel" 1 in the blue 
(red). We oriented the slit with PA = 105°, 133°, and 172° 
to go through the nuclei A and B, B and C, and A and C re- 
spectively, as labeled on the SDSS image shown in Figure Q] 
The total effective exposure time was 10800, 8100, and 8100 
s for A, B, and C, respectively. Standard stars G191B2B and 
HZ44 were observed for spectrophotometric calibration. We 
reduced the DIS data following standard IRAF 7 procedures 
dTodvill986h . We extracted one-dimensional spectra using 3" 
diameter apertures for the nuclei A, B, and C respectively. 
We applied a telluric correction from standard stars on the ex- 
tracted one-dimensional spectra. The median signal-to-noise 
ratio (S/N) achieved was 50-80 pixel" 1 . 

Figure [2] shows the resulting rest-frame one-dimensional 
spectra and our best-fit models for the stellar contin- 
uum and the continuum-subtracted emission lines, respec- 
tively. We measure redshift by fitting the continuum with 
galaxy templates produced by population synthesis models 

7 IRAF is distributed by the National Optical Astronomy Observatory, 
which is operated by the Association of Universities for Research in Astron- 
omy, Inc., under cooperative agreement with the National Science Founda- 
tion. 
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FIG. 2. — DIS one-dimensional spectra of the three nuclei of SDSS 
J1027+1749. We show the total spectra, the continuum-subtracted emission 
lines, and the continuum-and-emission-subtracted residuals offset by — 10~ 16 
erg s" 1 cirT 2 A "' . Plotted in red are our best-fit models for the continuum 
and emission lines, respectively. 



of iBruzual & Chariot! {2003 ) using the procedure described 
in detail in iLiu et alj (120091) . The templates have been con- 
volved with the stellar velocity dispersions er» measured over 
the spectral region of 4130-4600 A containing the G band 
A4304 A using the direct fitting algorithm of iGreene & Hoi 
( 2006) andlHo et alJ (l2009h following the procedure described 
in ILiu et all d2009l) . After subtracting the stellar continuum, 
we fit the emission lines simultaneously with Gaussian mod- 
els constrained to have the same redshift and width. 

3. RESULTS 

Tables Q] and [2] list our measurements of the host 
and emission-line properties of the three nuclei of SDSS 
J 1027+ 1749, respectively. The redshifts measured from stel- 
lar continuum fitting show that the components B and C are 
offset from A by ~ 450 and ^110 km s" 1 in velocity. Their 
stellar velocity dispersions er* are 200-300 km s" 1 , typical of 
massive galaxies at z ~ 0.1 of similar luminosities and sizes 
dBernardi et al.ll2003b . To securely separate the three compo- 
nents, we applied aperture photometry for them individually 
with 3" diameter (corresponding to 3.8 kpc) apertures using 
their SDSS images. They all fall in the red sequence on the 
u-r, M r color-magnitude diagram dBaldrv et al.ll2.004b . We 
have also performed surface brightne ss profile fitting using 
the GALFIT package (Pe ng et al.ll2010l) . using stars within the 
field to model the point-spread function (PSF). Assuming a 
five-component model for the entire galaxy (four Sersic pro- 
files for the three nuclei and the surrounding diffuse compo- 
nent, and an exponential disk for the tidal feature), we esti- 
mate their r-band Sersic indices as 1.4, 4.1, and 1.2 for A, 
B, and C, respectively. While the fit seems reasonable with a 
reduced \ 2 of 1.2, the results are highly uncertain given the 
limited resolution (r-band PSF FWHM ~ l."l) of the SDSS 
images and systematic uncertainties from model degeneracy. 
Hubble Space Telescope (HST) and/or ground-based imaging 
assisted with adaptive optics will help better determine their 
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TABLE 1 

Host-galaxy Measurements of the Three Nuclei of SDSS J1027+1749. 



r M r M, ;c u-r R e logM, logM,, c it* logM. 

ID SDSS Designation z (mag) (mag) (mag) (mag) (kpc) (M ) (Mq) (kms -1 ) (M Q ) 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) 

A J102700.40+174900.8 0.0668 16.73 -20.66 -22.4 2.48 1.1 10.4 11.2 265±6 8.6 

B J102700.55+174900.2 0.0664 17.00 -20.38 -23.1 2.46 1.6 10.1 11.2 193 ± 10 8.1 

C J102700.38+174902.6 0.0652 17.03 -20.30 -24.5 2.51 2.4 10.2 12.0 305±12 8.9 



NOTE. — Col. 1: ID of the three nuclei as labeled on Figure [T] Col. 2: SDSS names with J2000 coordinates given in the form of "hhmmss.ss+ddmmss.s"; 



Col. 3: redshift measured from stellar continuum fitting; Col. 4: SDSS r-band fiber magnitude; Col. 5: SDSS r-band fiber absolute magnitude; Col. 6: Col. 5 
after correction for internal dust extinction estimated using the Balmer decrement method; Col. 7: SDSS color from fiber magnitudes; Col. 8: effective radius 
from aperture photometry; Col. 9: stellar mass estimates from population synthesis modeling of the optical stellar continuum; Col. 10: Col. 9 after correction 
for internal dust extinction; Col. 1 1 : stellar velocity dispersion from model fits of the stellar continuum over the G-band A43 04 A region (Fi gureH). The quoted 
errors throughout this table are statistical uncertainties; Col. 12: black hole mass estimates assuming the M.-cr* relation of Tremaine et al. (2002) observed in 
local inactive galaxies. 



TABLE 2 

Emission-line Measurements of the Three Nuclei of SDSS J1027+1749. 



logL [0 iii] logL[oin]. c o-g as EW Ha log« e E(B-V) 

ID (L Q ) (L q ) (kms 4 ) [Om]/H/3 [Nn]/Ha [Oi]/Ha [Sn]/Ha Hct/H/3 (A) (cm" 3 ) (mag) 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) 



A 6.35 7.2 248±17 1.2±0.1 0.98±0.02 0.18±0.01 0.61 ±0.02 5.3±0.4 4.5±0.1 -2.0 0.61+g | 

B 6.38 7.7 211 ±6 0.52 ±0.03 0.66 ± 0.01 0.074 ±0.003 0.35 ±0.01 7.2 ±0.2 19.1 ±0.1 0.3+'| 0.93;* 

C 6.42 8.4 301 ±23 3.2±0.6 1.22±0.01 0.21 ±0.01 0.81 ±0.02 12.1 ±2.0 6.8±0.1 1.9^ 1.5^1 



NOTE. — Col. 1: same as Col. 1 in Table[T] Col. 2: [OIII] A5007 luminosity from Gaussian fits to the emission lines; Col. 3: Col. 2 after correction for 



internal dust extinction estimated using the Balmer decrement method as listed in Col. 12; Col. 4: gas velocity dispersion from Gaussian fits of the continuum- 
subtracted emission lines (Figure|2). The quoted 1 a errors throughout this table are statistical uncertainties; Cols. 5—9: intensity ratios from Gaussian fits of the 
emission lines; Col. 10: rest-frame Ho equivalent width; Col. 11: electron density estimated using the emission-line intensity ratio [S II] A6717/[S II] A6731. 
Col. 12: color excess estimated from the emission-line intensity ratio Ha/H/3 (Col. 9) using the Balmer decrement method, assuming the intrinsic case B values 
of 2.87 for T = 10 4 K )Osterbrock|[l989T) and the extinction curve of lCardelli ef"aTl f!^89l) with R v = 3. 1 . 

structural properties. We estimate stellar mass M» based on 
population synthesis modeling of the optical stellar contin- 
uum. The stellar mass ratio of the triple system A : B : C is 
2 : 1 : 1.3 before extinction correction and is 1 : 1 : 6 after. 
We estimate internal dust extinction from the emission-line 
intensity ratio Ha /H/3 using the Balmer decrement method 
(Osterbrock 1989) assuming the intrinsic case B values of 
2. 87 for T = 10 4 K (Oster brockl 19891) and the extinction curve 
o flCardelli et al.l (119891) with/?v = 3.1 . However, the estimates 
are highly uncertain due to the systematic uncertainty in the 
dust geometry. The three stellar comp onents are blended in 
Two Micron All Sky Survey (2MASS: ISkrutskie et al.ll2006l) . 
Higher resolution IR imaging is needed to better constrain 
their individual stellar masses. 

Figure [2] shows that all the three nuclei emit narrow lines 
with gas velocity dispersions of cr gas ~ 200-300 km s" 1 (Table 
O, consistent with <r* of their associated stellar components. 
Their emission line intensity ratios (Table [2) are all character- 
istics of AGNs or AGN-H n composites. The masses of the 
AGNs are estimated to be M. > 10 8 M© (Table [TJ from cr„ 
of their associa ted stellar componen t s 8 . Figure [3] shows the 
BPT diagram dBaldwin et alj 1 198 It IVeilleux &Osterbrockl 
11987b for the three nuclei. Based on [Nil] A6584/Ha and 
[O in] A5007/H/3, A is classified as a low-ionization nar- 
row emission region (LINER), B is a composite, and C 
is a Seyfert (or a LINER within uncertainties). A (C) is 
consistently classified as a LINER (Seyfert) based on diag- 

8 This still holds after correcting for the fact that M. is likely overestimated 
by ~ 0.2 dex us ing cr» measured in AGN mergers separated by a few kpc 
)Liuetal.ll2011J) . 
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FIG. 3. — Optical diagnostic emission-line ratios for the three nuclei mea- 
sured from the continuum-subtracted slit spectra (Figure |2}- Gray scales 
indicate number densities of 31,179 emission-line galaxies from the SDSS 
DR4 iKauffmann et al. 2003). T he solid curve is the e mpirical separation be- 
tween H n regions and AGNs iKauffmann et al. 2003), and the dashed curve 
is the theoretical "starburst limit" (Kewleyetal. 2006). Pure star-forming 
(SF) galaxies lie below the solid curve, AGN-dominated galaxies lie above 
the dashed curve, and AGN-H n composites lie in between. The dotted curve 
is the empirical separation between Seyferts and LINERs (Ho et al. 1997). 
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nostic ratios [S n] AA6717,6731/Ha or [Oil A 6300/Ha and 
O nil A5007/H/3 according to the criteria of iKewley et al.1 
2006). While LINERs and composites may be du e to stellar 
or sho ck heating rather than AGN excitation (e.g., iLutz et ail 
1999), X -ray obs ervation s show that ~ 80% of them do har- 
bor AGNs ("e.g.. lHoll2008l) . Because the intrinsic AGN lu- 
minosities of the three nuclei are obscured in the optical, we 
estimate the intrinsic hard X-ray luminosity ,2-iokeV using 
the [0 111] A5007 luminosity (corrected for internal extinction; 
Table [2j as a surrogate 9 , finding log(Lx,2-iokev/ergs _1 ) ~ 
42.4, 43.0, and 43.9 for A, B, and C, respectively. The 
non- detection of the s ystem in the ROSAT All-Sky Sur- 
vey dVoges et alJ 1 19991) is consistent with these estimates. 
The intrinsic [0 111] A5007 luminosities and Ha equivalent 
widths (Table [Hi are typical of z ~ 0.1 obscured AGNs (e.g., 
Kauffmann et al. 20031: iHol 120081) . The component B has a 
strong post-starburst component, as indicated by the strong 
Balmer absorption in our slit spectrum. This is also seen 
in its SDSS fiber spectrum which covers the higher order 
Balmer series and the Balmer break. B is also detected in 
the FIRST 1.4 GH z survey with an int egrated flux density of 
13.57 ±0.14 mJv dBecker et al.lll995l) . The galaxy is con- 
tained in the IRAS faint source catalog (Moshir et al.lll992b 
with a positional match centered on B, IRAS 10242+1804, 
with an IR (1-1000 ^m) luminosity of L IR ~ 1O 1L3 L esti- 
mated based on its flu x densities a t 60 an d 100 /im from IRAS 
using the formulae of lHelou etafl (11988b . 

4. DISCUSSION 

Based on dynamical friction timescale estimates 10 , the stel- 
lar components B and C will merge with A in t&ic ~ 40 and 
i~ 8 Myr, respectively, corrected for projection assuming ran- 
dom orientation for r p . The three SMBHs will inspiral with 
their individual stellar bulges at first, but will then decay to 
the center of the merged galaxy under dynamical friction with 
the stellar background in ~ 40-200 Myr 11 , given their mass 
estimates (Table [TJ. However, the actual evolution involving 
three comparable stellar components and the associated gas 
is likely more c omplicated. High r esolution kinematics ob- 
servations (e.g., lAmram et al.ll2007|) combin ed with tailored 
numerical modeling (e.g., iRenaud et al.ll2.01Qb are needed to 
unravel the merger history and future of this triple system. 
Regardless of the actual sequence of the mergers, the SMBHs 
may form a gravitationally interacting triple system, if the 
orbit of the third SMBH decays rapidly before the first two 

9 Measurements of Z,x,2-l0kev/i-[Om] f° r optica lly selected obscured 
AGNs range from a few to a few hundred (e.g., Heckman et a ll 120051; 
IPanessa et al. 2006). Here, we adopt the mean calibration of IPanessaeTal] 
(2001, log(Lx, 2 -iokev/ergs-') = 1.221og(L [0 , I i]/ergs- 1 )-7.34. 

10 Assuming singular isothermal spheres for the density profiles and circu- 
lar orbits, the dynamical friction timescale for a satellite galaxy with velocity 
dispersion ag inspiraling from radius r in a host galaxy with velocity disper- 

sion a H can be estimated as ffnc = fa ^ ( ^gprV ( 1 ^ £ ^) ' Gyr, 

where A is the Coulomb logarithm (Binney & Tremaine 1987). In A ~ 2 for 
equal mass mergers (Dubinski et al. 1999). We assume A as the host, and B 
and C as satellites, since A is at the apparent center of the galaxy. Assuming 
C as the host instead, which has the largest o* , yields similar estimates for 
'flic, ~ 20 Myr for A and C, and ~ 70 Myr for B and C, respectively. 

1 1 Assuming a singular isothermal sphere for the density distribution of the 
host galaxy, the dynamical friction timescale of a BH of mass M. on a circu- 
lar orbit of radius r can be estimated as ?f,i c = nr ( tt— I ,„„f H — r 1O ,. M0 

™ InA \5kpc / 200 kms -1 M. 

Gyr, where In A ~ 6 for typical values (Binney & Tremaine 1987). We es- 
timate that ffnc ~ 200 and 40 Myr for the SMBHs in B and C, to reach the 
center of the merger products of B and C with A, respectively. 



SMBHs merge. 

The fraction of kpc-scale pairs of AGNs is /kpc, double ^ 
3 x 10~ 3 corrected for SDSS spectroscopic incompleteness 
based on the number of A GN pairs with r n < 10 kp c found in 
our p aren t AGN sample (|Liu et al.ll201 lbt see also iLiu et al.l 
120101 and lShen et al.ll201 ll for a similar estimate of /k pc , double 
from a complementary approach to identify kpc-scale AGN 
pairs, based on the selection of AGNs with double-peaked 
narrow emission lines). This is a lower limit due to the lim- 
ited imaging resolution of SDSS. Out of the seven kpc-scale 
triple AGN candidates, we have obtained slit spectra for five 
systems, including SDSS J1027+1749. The nature of the 
other four systems seems less clear, due to alternative scenar- 
ios involving star clusters or tidal knots, the results of which 
will be presented in a future paper. Our result suggests that 
the frequency of kpc-scale triples in optically selected AGNs 
at z ~ 0.1 is Apcripie > 3 x 10" 3 x ^ x \ « 5 x 10" 5 . If 
we assume that AGNs in successive mergers are independent 
events with short duty cycles, the expected frequency of kpc- 
scale triple AGNs based on the observed frequency of dou- 
bles would be ~ /k pc douMe > 10~ 5 . Our observed frequency of 
triples implies a higher probability than this, which may sug- 
gest that AGNs in mergers are somewhat correlated, although 
better statistics is needed to draw any firm conclusion. 

While our slit spectra of the three nuclei suggest that they 
all host AGNs, optical identification alone is inconclusive. Al- 
ternatively, there could be only two or even just one active 
SMBH, ionizing all the three gas components in the merging 
system. To discriminate between these scenarios, arguments 
based on the ionization parameter and effe ctive size of the 
narrow-line regio n are not very useful (e.g., Liu et al.l 120101 : 
IShen et alJl201lh . bee ause the separation between the nuclei 
are not much larger than the individual size of the gas-emitting 
regions, and the electron density measurements are highly un- 
certain (Tabled. Chandra X-ray observations could help pin 
down the triple AGN nature of the system. 
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